The Rashba Effect within the Coherent Scattering Formalism 
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The influence of spin-orbit coupling in two-dimensional systems is investigated within the frame- 
work of the Landauer-Biittiker coherent scattering formalism. This formalism usually deals with 
spin-independent stationary states and, therefore, it does not include a spin-orbit contribution to 
the current. In this article, we will rederive the coherent scattering formalism, accounting for the 
changes brought about by the spin-orbit coupling. After a short review of the features of spin- 
orbit coupling in two-dimensional electron gases, we define the creation/annihilaton operators in 
the stationary states of the spin-orbit coupling Hamiltonian and use them to calculate the current 
operator within the Landauer-Biittiker formalism. The current is expressed as it is in the standard 
spin-independent case, but with the spin label replaced by a new label which we call the spin-orbit 
coupling label. The spin-orbit coupling effects can then be represented in a scattering matrix which 
relates the spin-orbit coupling stationary states in different leads. This scattering matrix is calcu- 
lated in the case of a four-port beamsplitter, and it is shown to mix states with different spin-orbit 
coupling labels in a manner that depends on the angle between the leads. The formalism is then used 
to calculate the effect of spin-orbit coupling on the current and noise in two examples of electron 
collision. 
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I. INTRODUCTION 



Coherent electron transport through nanostructures in 
cryogenic two-dimensional electron gas systems is inher- 
ently a quantum mechanical phenomenon. Several exper- 
iments have demonstrated certain aspects of this quan- 
tum behavior, whether it relies directly on the wave na- 
ture of the electron and can be probed through a cur- 
rent or conductance measurement {e.g., quantized resis- 
tance in the quantum Sail effectu, conduction modes of a 
quantum point contactfiu), or on the particle nature and 
quantum statistics of the electrons and can be probed 
through a noise measurement (e.g., Hanbucy Brown and 
Twiss-type experimented, electron collisionu, observation 
of the.jfractional charge in the fractional quantum Hall 
effectOn). In addition, 2D electron gases could be used to 
study the fundamental non-local features nf rui ant-uia me- 
chanics through electron entanglementLtE30t3'LJO. All 
these experiments can be successfully explained within 
the coherent scattering formalismli^O, a theoretical tool 
describing coherent and non-interacting electron trans- 
port. This formalism relies on spin- independent station- 
ary states in the leads of the device and, therefore, de- 
scribes spin-independent transport. Although it is pos- 
sible to add a local spin-dependent effect directly in the 
scattering matrix, it is not possible in general to take 
into account spin effects occurring over the whole sys- 
tem. One potentially important spin effect occurring in 
the leads of the conductor is spin-orbit coupling. Any 
electric field in the reference frame of the laboratory gen- 
erates a magnetic field in the moving electron reference 
frame, coupling the electron's orbital degrees of freedom 
with its spin. One can find several sources of electric 
fields in semiconductors. In three dimensional crystals, 
the periodic crystal potential generates the Dresselhaus 



cffectS, which induces a spin-splitting of the conduc- 
tion band that is proportional to fc 3 . In two-dimensional 
systems, the dominant term results from the asymme- 
try of the ia-plane confining potential. First introduced 
by RashbaLS, this effect causes a spin-splittingEj propor- 
tional to k, and it depends on the strength of the applied 
electric field. Recently, there has been a growing interest 
in electronic devices which rely on the spin properties of 
the electrons. These spin-dependent devices may be in- 
fluenced by spin-orbit coupling effects, or may even rely 
on it as, for example, in a coherent version of a spin- 
polarized field effect transistorc£lEJ. Therefore, it may 
be useful to include spin-orbit coupling in the coherent 
scattering formalism. 



II. SPIN-ORBIT COUPLING IN 
TWO-DIMENSIONAL ELECTRON GASES 

The spin-orbit coupling is governed by the following 
spin-orbit HamiltonianE3, which is obtained using an ex- 
pansion in v I c of the Dirac equation 



H 



so 



(2to c) ; 



-W(<rxp), 



(1) 



where too is the free electron mass, p is the momentum 
operator, <r = (d x , cf y ,cf z ) are the Pauli spin matrices , V 
is the electrostatic potential, and V is the gradient oper- 
ator so that — W is the electric field. Electron transport 
in the presence of an electric field results in a spin-orbit 
effect which couples the electron spin and orbital degrees 
of freedom through the <x x p term. Here, we will ne- 
glect the bulk effect arising from the periodic crystal po- 
tential (Dresselhaus effect) and consider only the effect 
caused by the asymmetry of the confining quantum well 
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(Rashba effect), as it is stronger than the Dresaelhaus ef- 
fect in most two-dimensional heterostructuresE3oc3. Let 
z be the direction of confinement, perpendicular to the 
plane of motion. The asymmetry of the confining po- 
tential along the z-direction results in a non-zero electric 
field along the z-axis (E = — Eq u z throughout the arti- 
cle, where u z represents a unit vector in the z-direction). 
The spin-orbit coupling term of the Hamiltonian can be 
written 

a . A A . , „ d „ d . , . 

Hso = j\p X P)z = ioi{a y - g x 7^— ) , (2) 

x y 

where a is the spin-orbit coupling constant and de- 
pends on the strength of the electric field. It takes 
values in the range 1 to 10 x 10~i2, eV cm for,— a 
large variety of systemSj-ZJoAs/GaSbcJ, InAs/AlSbeJ, 
In^Gai-^As/In^Ali-^AsE^EZl and GaAs/Al x Gai_ x AsEa) 
depending on the shape of the confining well. Using the 
standard effective mass approximation, we can deduce 
the system Hamiltonian as the free Hamiltonian plus the 
spin-orbit coupling Hamiltonian 



H = 



Px 



Py 



2m 



-(VyPx - <7xPy), 



(3) 



where m is always taken to be the effective mass. Since 
the operators p x and p y commute with H, we can search 
for eigenstates of the form 



(4) 



where |f) and ||) label the up and down states of the 
z-component of the spin. We can now diagonalize the 
Hamiltonian 



H 



h A k A 

2m 



ak v — iak x 



ak y + iak x 
h 2 k 2 

2m 



(5) 



in this spin subspace. The eigenvalues are E(k) — ± 
ak, and the associated eigenfunctions are 

p i(k x x+k y y) 

= [ e * 3 It) + e-*-* ||)] (6) 



V2 

p i(k x x+k y y) 



IT) 



14)] (7) 



where 9 is the angle between k = (k x , k y ) and the y-axis 
(see Fig. Q). The electrons feel a virtual magnetic field 
in the 2D plane in a direction perpendicular to k. The 
spin is aligned or antialigned to this field (see Fig. [I]) , 
so that 0(k,Si Wj5 v) = vr/2, and 0(k,Si Wj5 p = -7r/2, 

where 0(k, Si^_\) is the angle between k and the spin 
S in the state I )• The amplitude of the magnetic 
field depends on the velocity of the electron and van- 
ishes for k = 0, preventing a possible spin polarization 
in the system. The spin-orbit splitting is usually small 
compared to the kinetic energy of the electrons ( 0.15 
meV < akf < 1.5 meV for Ef — 14 meV). Following 
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FIG. 1: Direction of the virtual magnetic field and spin. An 
electron with wave vector k in the presence of an electric field 
E = —Eq u z feels a virtual field B(k, E) perpendicular to k. 



Rcf. [l5|, we will introduce a transverse confinement in the 
leads of the conductor, allowing us to address the longi- 
tudinal transport modes for each given transverse mode. 
Focusing on one lead, 7, we will make two simplifying as- 
sumptions. First, we consider only a single independent 
transverse mode. Second, we neglect the ID SO cou- 
pling effect that this transverse confining potential could 
create, since, to our knowledge, there is no experimen- 
tal verification of this effect, and iLis estimated to be 
much smaller than the Rashba effectcj. Within these ap- 
proximations, we can use our previous analysis to deduce 
the eigenstates and the associated energy dispersion di- 
agram (see Fig. ^|), with k lying in the direction x 1 of 
the lead (making the angle 1 with the y axis). We will 
now introduce three labels for the eigenstates that will 
prove useful in writing the creation and annihilation op- 
erators of these states, e = a or b labels the direction of 
propagation from the sign of the group velocity v g (a if 
sgn(v g ) > 0, b otherwise). The parameter k labels the 
longitudinal mode wavevector. The SO coupling label 
(7 = ± designates the two different branches of the en- 
ergy dispersion diagram for a given k (+ branch and — 
branch in Fig. ||). Using these labels, we find the follow- 



|b, k 




[71 + states k = 
[71 - states e„ = 



ma 

2h 2 



FIG. 2: Energy dispersion diagram with SO coupling. The 
—rjako < k' < 00 range (a states) is deduced from Eq. (^) 
with k' — k. The —oo<k'<r]ako range (b states) is then 
deduced from mirror symmetry (k' — —k). 
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ing eigenstates and eigenvalues of the system from Eq. 
and® 



III. SPIN-ORBIT COUPLING AND COHERENT 
SCATTERING FORMALISM 



\e,k,a) = $(y 7 )__[ e ^|t)+ W^ll)] 
hk 2 

E = — h rj a ak — f] a k < k < oo (8) 

2m 

where rj c = +1(— 1) when e = a(b), and r\ a — +1(— 1) 
when a = +(— ). &(y-y) is the normalized transverse 
wavefunction for the transverse mode under consider- 
ation. By convention, k is only taken in the range 
—V&ko < k < oo for both the a and b states. We 
use r] e = ±1 to parameterize explicitly the appropriate 
wavevector range for the e = a(b) propagation direction 
in the eigenstates' phase factor e iT, ' kx . The eigenvalues in 
Fig. U can be deduced by mirror symmetry. This conven- 
tion allows us to track the propagation direction through- 
out the calculation. We also notice in Fig. || that, for 
a given energy, the corresponding + and — states with 
same label e do not have the same wavevector: 

ATT} cy 

k{E,-)-k(E,+) = Ak = 2ko = ^. (9) 

h 

The states with a = ± have their spin perpendicular to 
the direction of propagation v g and in opposite direc- 
tions, so that 0(v g ,S^.-)) — k/2 and 0(v g ,S|fc i+ )) = 
— 7r/2 (see Fig. ph. It is important to notice that the 
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FIG. 3: Direction of propagation and direction of the spin 



a = +(— ) states do not coincide with the | V E , £; + ) 
states, as it is v g , and not k, which determines the di- 
rection of the spin (e.g., in Fig. the spin of the + state 
is perpendicular to and on the right side of v g ). For ex- 
ample, in Fig. ||, the two definitions are consistent for 
states with positive (negative) k! and group velocity v g . 
However, they are inconsistent for positive (negative) k' 
and negative (positive) v g . We choose here to link the 
direction of the spin to v g (rather than k') through the 
label a = ± (rather than E±), because it is the group 
velocity which determines the direction of the current. 
Note that in the case without SO coupling, no differen- 
tiation between v g and k! is necessary since they always 
share the same sign. 



The Landauer-Biittiker formalism relies on a second 
quantization formulation of quantum mechanics. The 
current is expressed as a function of the field operator, 
which is expanded in the basis of the scattering states 
using the operators that create and destroy electrons in 
the leads of the conductor. A scattering state is a coher- 
ent sum of an incident wave in one lead and the outgoing 
waves it generates in all the other leads. The amplitude 
of each outgoing wave is given by the scattering matrix 
S, whose elements depend on the properties of the scat- 
terer. As our final goal is to study the fluctuations of 
the current, it is important to find a simple way to ex- 
press the time dependence of the current operator. This 
leads us to choose the stationary states in each lead for 
the purpose of expanding the field operator. In our red- 
erivation of the coherent scattering formalism, we will 
follow the same analysis as in Ref. [l5|, but introducing 
two important changes: we will rederive the expression 
of the current operator as a function of the field operator, 
adding a new contribution coming from spin-orbit cou- 
pling, and we will expand the field operator in the new 
stationary basis of the SO coupling Hamiltonian. 



A. Calculation of the current operator 

One way to find the current operator for non- 
interacting electrons is to use the single-particle 
Schrodinger equation 



0* 



9* 



dt 



2m dx 2 dy 2 



dx 



dy 



( 10 ) 

where VP is a two-components spinor. The current oper- 
ator j can be extracted from the conservation of charge 
equation ^[e* 1 ^*] + V.j = 0; using Eq. @ and its 
adjoint, 



eh r +5* 
; - 



Jx = 



Jy 



2mi 
eh 
2mi 



dx dx 







ea 

~h 

ea 



^CTj,* (11) 
*V S * (12) 



dy dy 

We can identify the usual kinetic term of the current 
density 

eh 



JK 



2mi 



(13) 



Note that SO coupling adds a new contribution propor- 
tional to a that we call the SO coupling current density 
jso, 



Jso 



en 



tfVstfUy. (14) 

h 



In the framework of second quantization, ^ becomes a 
field operator, and we have to expand it using a conve- 
nient basis: the stationary states. 
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B. Expansion of the field operator 

The first step is to define the creation (annihilation) 
operators for the SO coupling stationary states. We first 
introduce the creation operators in the old spin basis: 



a 



7feT $vki) crea tes an incoming (outgoing) electron with 
spin up (down) in the first mode of lead 7 with momen- 
tum k, in the state |o, k, f) (\b, k, J.) ), where |e, k, T) = 
e" ukx '< \\). These standard operators satisfy the anticom- 
mutation relation a 7i fc iS , k , s ,]_|_ = Skk'S-ypfiss' , where 
s labels the spin of the electron. Through Eq. (|J) we in- 
troduce , which creates an incoming electron in lead 
7 with momentum k and spin-orbit coupling label a (not 
to be confused with the spin s) satisfying the relation 



1 



V2 



(15) 



Knowing the anticommutation relationships in the spin 
basis, we can calculate it in our new stationary basis to 
be 



Ok. 



1," > ^!',/3,<t 



= Skk'SafsScrcr' , 



(16) 



where we used the relation rj a rj a i = 2S aa i — 1. This result 
justifies our labeling of the SO coupling eigenstates and 
the use of the spin-orbit coupling label a instead of the 
spin. The stationary states form a complete basis, and 
we can use them to expand the field operator in lead 7. 
From now on, we will choose x 7 and y 7 as the reference 



axis (so that 6L 



specifically tracking all of the 



angular dependence in the scattering matrix. 



(T,e k— — 'r] cr k a 



(17) 



where e = a(b) is the annihilation operator in the incom- 
ing(outgoing) states and x% is a two-component spinor 
obtained from Eq. (||), setting # 7 = ^ 



e ir, c kx 

7z 



l+i 
1. 



(18) 



As shown in Fig. its spin depends on the SO coupling 
state (a) and on the direction of propagation (a or 6) . At 
this point, it is more convenient to go from a discrete sum 



in £;-space to a continuous integral in energy by defining 
the operators (E) , which create electrons in the energy 
continuum. As our labeling defines a one-to-one corre- 
spondence between (e,k,a) and (e, E, er), we can define 
unambiguously a 7(T (i?) = \J p{E)d\ fk(j} so that 



a^ a {E) , a^A E ')} =6^6 aa ,d(E-E') (19) 



where p(E) is the density of states at energy E . A priori, 
p(E) could depend on the spin-orbit coupling label ct, but 



(°xy) 



FIG. 4: Direction of the spin of \% 

one can check in Appendix [a] that it is not the case. Re- 
placing the discrete sum over k by an integral over E and 
introducing the new operators and their time-dependence 
a(E,t) = e'^^E), we find 



*(x,?/,i) 7 



-y 



30 jp 

Eg ^hVg(E) 



e i Ve k(E.*)x e ~i§t (20) 



with Eq = — %jj£r, and k has been replaced by k(E,a) 
to remind us that, for a given energy, k depends on the 
spin-orbit coupling label a. 



C. Current in the spin-orbit coupling basis 

The current in lead 7 is given by 

I a = I dyj x = / dy(j xK + jxso) = Ik + Iso- (21) 



We will begin by calculating the usual kinetic term Iff. 
Using Eq. (O), we have 



= EE / dEdE* i K i a ,(E,E>) i+ a (E) e\ a ,(E') e*^. (22) 



00 e e 



I K iAE,E') = l X°X°> [ Ve k(E,a)+ Ve ,k(E\a')} j^W)-***,*))*. ( 23 ) 

2mi hy/v g {E)vg{E') 
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We notice from Eq. ( p2| ) that the u> frequency of the cur- 
rent is given by E — E'+hcu. Following Ref. [15| we calcu- 
late the current and noise in the zero-frequency limit, and 
make the approximation k(E' , a) « k(E, a) (as E s» E') 
for electrons having the same SO coupling label, and, 
using Eq. (||), k(E',a') = k(E,a) + i] a Ak for electrons 
having different SO coupling labels. From Eq. (|l8|), we 
also have, (with S e ^ e ' = 1 if e ^ e', otherwise) 

XctXct' = <W<W + ^'^ct'- (24) 



This result is consistent with the fact that only electrons 
with the same direction of propagation and the same 
spin-orbit coupling label, or opposite direction of prop- 
agation and opposite spin-orbit coupling label, have the 
same spin. 



In the a = a' case, implying e = e' and k(E' , a') ss k(E, u), we find 



7 if ctct'IA-^ ) l<W - 1eT 



7i 



■ma 



(25) 



In the a ^ a 1 case, implying e 7^ e' and k{E',a') w er) + r^^, so that r) e k(E,a) +r^k / (E,a') = -r) e rj a ^, 
we find 



Vo 



h k(E,+) + f£ 



-ir] € (2k(E,<r)+r] tT Ak)x 



(26) 



If we only consider the kinetic term of the current, a non- vanishing contribution for e ^ e' leads to terms like and 
wd in the expression of the current, corresponding to electrons propagating in opposite directions. 

Now, let us study the contribution of the spin-orbit coupling current Iso- Using Eq. ( jl4| ) to calculate I£o ao 1 m 
Eq. (p2|) with k^so, we find 



ISO aa'( E i E ') =~T 



h ^v g (E)v g {E>) 



Act °yXa' e 



(27) 



Using ije'V = ??c?7ct, we have from Eq. (|l|) ct^Xct' = -WX^/ and Xct °2/Xct' = (<W <W + 6e?e> <Wct')> so 

that 



h k(E,+) + ^ 



A , A . _i_ p -iVe(2k(E,cr)+r] a Ak)x X , , X , 
u ee' u ctct' 1 e u e^e' u ct^ct' 



(28) 



We can now calculate the value of the total current from 
Eq. ©, @ and @ 

1 (25, E') = I£ i a , (E, E') + IiS'^ (E, E') 

= Ve^ 5 ee > <W- (29) 

Although the spin-orbit coupling mixes electrons having 
different directions of propagation when we consider only 
the kinetic term of the current, these contributions cancel 
when we add the spin-orbit coupling current Iso , and we 
find the standard formula for the current operator, 

I = lYl f dEdE' [a\ a {E)a ia {E') 

-b\ a (E) b ia (E') } e^K (30) 

In this formula, the definition of the d(b) states as states 
with positive (negative) group velocity, and not necessar- 



ily positive(negative) wavevector, is consistent with the 
fact that they carry the current in opposite directions. 
More importantly, the final expression of the current is 
similar to the one found in Ref. but with the spin in- 
dex replaced by the spin-orbit coupling label a — ±. The 
spin related to this new index depends on the direction of 
propagation, that is, the angle of the lead 9 1 . Therefore, 
we expect a Q 1 dependence in the scattering matrix relat- 
ing the spin-orbit coupling states in leads with different 
orientations. As an example, we investigate the scatter- 
ing matrix in the case of a four-port beamsplitter (two 
input leads, two output leads) used in electron collisions. 
The scattering matrix relates the outgoing states in the 
outputs (b states) to the incident states at the input (a 
states) 

f3a' 
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In the spin-independent transport the beamsplit- 

ter does not act on the spin degrees of freedom, the scat- 
tering matrix does not mix different spins (the scattering 
matrix is diagonal in each subspace of spin: spin-up and 
spin-down). In this case the reflection and transmission 
coefficients do not depend on the spin. However, when 
we include spin-orbit coupling, the situation is more com- 
plicated, because the spin associated with the spin-orbit 
coupling label a is different in each lead (the leads having 
different orientations). The conservation of the spin at 
the beamsplitter then implies that we have a mixing of 
the spin-orbit states (off-diagonal elements in the scat- 
tering matrix) , and this mixing becomes more important 
when the angle between the leads increases. A beam- 



FIG. 5: Four-port beamsplitter 
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splitter, with all leads oriented at the same angle 8 (the 
case illustrated in Fig. ||) is investigated in Appendix 
and the following scattering matrix is found after some 
simplifying assumptions 

















\ 



i r sin ( 
t 




i r sin ( 
r cos 8 



t 



We note that the unitarity of the S matrix requires 
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r cos 8 
-i r sin( 




t 

-i r sin( 
r cos 8 







) 


ai- 




{ ° 2+ J 




\ 0,2- I 



s 



M 2 + |t| 2 



(32) 



R + T = 1 and rt* + tr* = 0. 



IV. SPIN-ORBIT COUPLING AND NOISE IN 
ELECTRON COLLISIONS 

We will now use the expression of the current oper- 
ator and the scattering matrix to calculate the current 
and noise in two experiments using the four-port beam- 
splitter setup. The first is an unpolarized electron colli- 
sion, which has already been realized experimentally (see 
Ref. ||). The second is a spin-polarized electron collision. 

A. Unpolarized electron collision 

This experiment is a collision at a 50/50 beamsplit- 
ter of electrons from leads I and 2 which are biased at 
Ef + eV. At zero temperature and without spin-orbit cou- 
pling, each electron with a specific spin in lead 1 has an 



identical partner in lead 2. The Pauli exclusion principle 
then forbids us to find these two electrons in the same 
output. We expect a total suppression of the partition 
noise. With spin-orbit coupling, we do not expect this 
result to be modified. It relies only on the full occupa- 
tion of the energy eigenstates, which is not modified at 
zero temperature by SO coupling. We will confirm this 
by studying the collision of four electrons at the same 
energy. The initial state is: 

1*0 - 4>,-,x( E ) 4,+, 2 (£) 4<,-. 2 (£) io> , 

(33) 

From now on, we will drop the E, k and k' labels for the 
energy and the momentum of the electrons, associating a 
+ (— ) SO coupling state with k(k'). The output state is 
found using the new scattering matrix with SO coupling 
using Eq. Q and Eq. @ 
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r [ cos 9 St_ 3 + % sin ( 



P , 1 + 1 si 



r [ i sin 9 6^ 3 + cos 9 6l 3 



t b\ 3 + r 



[ cos 9 — i sin 9 61 4 ] 



tbl 



[-i sin 9 6^ 4 + cos 6» &1 4 ] |0) 



= bl 3 ll 3 blJU\o) 



(34) 



where we have used the relation r 4 + t A — 2 r 2 t 2 = 1. 
As expected, a full occupation of the two input states 
at energy E leads to a full occupation of the two output 
states as well. A calculation of the noise for the fully 
occupied outputs would show a complete suppression of 
the partition noise. The only difference with the spin- 
independent transport case is that the colliding electrons 
at the same energy do not have the same momentum (k 
k'). Therefore, the non-ideality in the noise suppression 
observed in Ref. ^ cannot be caused by the SO coupling. 



B. Spin-polarized electron collision 

In the previous example, an unpolarized collision can- 
not show any SO coupling effect, because, starting with 
two quiet sources of electrons with fully occupied en- 
ergy levels, all the transport statistics are governed by 
the Pauli exclusion principle independent of the SO cou- 
pling. However, this is not the case for a spin-polarized 
collision. As the standard spin basis is not the stationary 
basis, and as the SO coupling scattering matrix mixes dif- 
ferent SO coupling states, some of the partition noise is 
recovered, depending on the SO coupling constant a and 
on the angle between the colliding leads. We will con- 



firm this by calculating the non-equilibrium noise (of the 
electrons above the Fermi sea) between Ef and Ef + eV, 
assuming that the process of polarization does not effect 
the transport in the conductor (for example, there is no 
magnetic field in the conductor) . The input state is then 
made of all the electrons with the same spin, for example 
spin-up, between kf and fc max = k(Ef + eV). Follow- 
ing Ref. [l^, we present the initial state and the current 
operator using the following notations 



I*,;) = 



I a (t) = 



AOS 



f[ a f 1T (fc)4 T (fc)|0) (35) 

k=kj 



\cts 
j3~yaa' 



s aa' j3*y 

x a^(£)<V(£') (36) 
S*£ S S™, (37) 



Trkt 



We then 



where k ma x is given by Ef + eV 

use the current fluctuation 6I a (t) = I a (t) — ( I a (t) ) to 
calculate the fluctuation correlation function. 



(si a (t)SI a (0)) = e -Y, E E JdEdE'dE''dE'"e tlS ^A^ aal A^ a ,, a ,,, A^f "'" (38) 

ss' aa' a" a'" {3-y5£ 

a^ c ' v " ^ ($ a (E)s h A&)$A&')kM&")) - (aj a (s)v(^')) (&tAE")k*">(E'")) (39) 



Although we derived an expression for the current oper- with D a (k) = dE j k k ' a ' > = + rj a a. Defining II 

ator in the SO coupling basis, it is mathematically more ^jD G (k)D a , {k')D a „ (k")D a ,„ (k'"), we have, 
convenient in this case to express all the operators in the 
standard spin basis 



^2D a {k) 

(40) 
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. (7(7 (7 a 



1 

in 



5p(5 y s S(k E ,a - k E iii,a>") S(k E \cr> - k E " ,<r") n^(kE,a) [1 - n 7 i(k E ' ,„')] 

■■■ + r}<T'r}<T»Sl3(Sy6 5{k E , a - kE>" ,cr"i)?>{kE> ,a' - k E " ,a") n f ^(k E ,a) 



(41) 



where k E ,u = k E ' yCr < for E = E' + ak(r] a > — i] a ), so that 5(k EyCr — k E '",a'") — D cr (k)S[E — E'" + ak{r\ a — ?7 CT /»)]. After 
integration over E and E and making the approximation D+(k) ss -D_(/c), we find 



dE"dE"'A™'£'°" 



5pC,5jS n/3l(k E} a) [1 - n 7 l(kE',a-')] + Tla'r\a" S 7 s npi(kE,a) 

I 



We can replace D-{k) by D+(k), because 

gso(fc/) 



1 



1, which follows from S-r < 



D-{k) = 

Qfej: _ 



<C 1 as the SO coupling is small. We will now 
calculate the power spectral density at zero frequency in 
output lead number 3, £33(0), 



S33M 



dt e iwt (SI 3 (t) SI 3 (Q)) . (43) 



From Eq. ( J38| ) and Eq. (|42|), we find 

„2 



3s' 



• Sii::l)i 27; X! E X ^Iw • 

ss' (7(7' a" a'" ,3,7—1,2 



x dE npi(kE,o)[l - n~(]{k E ,a>) + f} a >f} a »] (44) 



(42) 



where 7 = 1 or 2, since we ignored scattering from output 
to output. np^(k E ,tr) is the number of electrons in lead 
P with spin up and momentum k E>a 

npt{k E , a ) = 

if E < E f + r] a ak f 

1 if Ef + J] a ak s < E < E f + eV + r\ G ak max . 
if Ef + eV + r] a ak max < E 

Assuming that the SO coupling is weaker than the bias 
voltage, that is, Ef + akf < Ef + eV — ak max , we have 
np r (k E -)[l - n 7T (fcE )+ )] ^ for E s - akf < E < Ef + 
akf and n^{k E .+)[l — n^(kE,-)] ^ for Ef + eV — 
ak max < E < E f + eV - ak max . Therefore, Eq. @) 
becomes 



s 33 (o) = |i E E 4W' [ eV + 5 °+ 5 °'- 2 ak ^* + 6 «- s °'+ 2 ak A- ( 45 ) 

ss' aa' a" a'" /3,7 



Given the scattering matrix calculated in Eq. (|32j), we 
find 

E E E^V A lv"VV=0, (46) 

ss' aa' a" a'" {3,j 

and 

E E E4W' A% a ,, <T ,,,5 rT+ 8 rT ^=ATRsvc?e, 

ss' aa' a" a'" {3,*y 

(47) 

so that 

S 33 (0) = \ i^kf + ak max ) T(l - T) sin 2 0. (48) 

After some calculation, we find the current and Fano fac- 
tor to be 



e 2 V 



(49) 



F 



S 33 (0) 
(h) 



4eT(l-T) sin 2 



q(fc/ + Awe) 



(50) 



This result reveals two typical features of SO coupling: 
first, the noise is proportional to the SO coupling con- 
stant, and second, it depends on the angle between the 
input leads 1 and 2. One can check that the obtained 
noise is identical to the partition noise for two indepen- 



dent leads where only a fraction 



1 a{kj+k max ) - 
eV 



of the elec- 



trons is colliding, and for which the noise is modified by 
the factor sin 2 9. This result can also be explained using 
Fig. [| In the inputs, the spin-up electron states between 
kf and k max are fully occupied. However, in terms of 
the SO coupling states, only the states between E(kf, +) 
and E(k max , — ) are jointly occupied, giving no contri- 
bution to the noise as in the unpolarized case. Between 
E(kf, — ) and E(kf, +), only the — states are filled, and 
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k 



without SO coupling with SO coupling 

no splitting splitting 

|^>j + states occupied only 
[>Bj + and - states occupied 
- states occupied only 



FIG. 6: Energy diagrams with and without SO coupling 



between E(k max ,—) and E{k m ax,+), only the + states 
arc filled. These states do not have the same spin in leads 
1 and 2 (their overlap is cos 2 9) , so that the Pauli exclu- 
sion principle does not fully suppress the noise, and the 
classical partition noise is partially recovered through the 
factor sin 2 = 1 — cos 2 9. In this treatment, we have as- 
sumed that the SO coupling is small (tZitLd^mssl < i) ; so 
that E(kf,+) < E(k max , — ). But, a typical value of the 
applied bias voltage is eV = 0.1 meV, which is smaller 
than the typical values of energy splitting caused by the 
Rashba effect (0.3 meV < 2akf < 3 meV). In this case, 
all the electrons contribute to the partition noise, and we 
have F = 4 eT(l — T) sin 2 9. Depending on the angle 
between the leads, we can go from full noise suppression 
to the classical limit of the partition noise. The modifi- 
cation of the noise caused by SO coupling in this collision 
experiment can in principle be measured, even for small 
values of the SO coupling constant. 



replaced by the SO coupling label, indicating the align- 
ment of the spin on the virtual magnetic field caused by 
the SO coupling. The main differences introduced by the 
SO coupling then arise in the calculation of the scatter- 
ing matrix relating the stationary states of SO coupling 
in different leads with different orientations. The direc- 
tion of the virtual magnetic field (and, consequently, the 
direction of the spin) depends on the direction of propa- 
gation which is different for each lead in general. There- 
fore, the scattering matrix is shown to mix states with 
different SO coupling labels, and the strength of this mix- 
ing depends on the angle between the leads. The effect of 
SO coupling on the current noise was then investigated 
in two examples of electron collision. In the unpolar- 
ized electron collision example, it is shown that the SO 
coupling does not modify the noise; this case is entirely 
determined by the Pauli exclusion principle. In contrast, 
the polarized case exhibits a contribution to the noise 
caused by SO coupling, which is proportional to the SO 
coupling constant a and depends on the angle between 
the leads. A polarized electron collision experiment pro- 
vides, in principle, another way to measure the strength 
of the Rashba splitting energy. Finally, this new formu- 
lation of the current operator can be applied to other 
coherent scattering experiments in which one wants to 
investigate or incorporate the effects of SO coupling. 
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V. CONCLUSION 

We have studied the influence of SO coupling in the 
framework of the Landauer-Buttiker formalism. A short 
review of the effect of SO coupling in 2DEGs (and more 
precisely of the Rashba effect) has reminded us of its main 
features: spin-splitting proportional to k, and stationary 
states of the spin perpendicular to the direction of prop- 
agation. The electron feels a magnetic field perpendic- 
ular to the direction of propagation, with an amplitude 
proportional to the velocity. We have then included the 
effects of SO coupling in the Landauer-Buttiker coherent 
scattering formalism. The SO coupling gave rise to two 
important modifications. First, the addition of a SO cou- 
pling term in the Hamiltonian modifies the expression of 
the current operator, and an extra term directly related 
to SO coupling has to be included. Second, the expansion 
of the current operator is in the basis of the stationary 
states of SO coupling. The final formula for the current 
operator was found to be identical to the one derived in 
the spin-independent transport case, but with the spin 



APPENDIX A: DENSITY OF STATES WITH 
SPIN-ORBIT COUPLING 

In this appendix, we calculate the density of states and 
show that it does not depend on the SO coupling label 
<7, as is suggested by the symmetry between the + and 
— states in the energy dispersion diagram (see Fig. ^), 



h 2 k 2 (E,<r) 

E = ±-L-L +r}<r ak{E,o) 

2m 



(Al) 



\dE{k,a) hk(E,a) a . . „. 



From Eq. (O) we deduce 



ma 



k{E,a) = k(E,+) + (r ](T -l)- T (A3) 

n 

, „s hk(E, +) a . . , 

v ga (E) = y m ' +- = v g (E) (A4) 



P(E) 



1 



hv g (E) 2irhv g (E) 



(A5) 
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which is independent of a. 



APPENDIX B: SCATTERING MATRIX IN THE 
FOUR-PORT BEAMSPLITTER CASE 

Here, we determine the scattering matrix for the four- 
port beamsplitter described in Fig. ^. The beamsplitter 
is very simply approximated by a potential barrier at V = 
Vq of length L. The plane is then divided into three areas 
of different potential (see Fig. ^) in which the solution 
of the Schrodinger equation is known. Starting with an 



I k r (E,+), +> 



Reflected wave 

\ (lead 3) 


Beam-Splitter : 
Potential barrier at V = V 

Refracted wave 




ru 




Transmitted wave 

(lead 4) 

<°x> 

► 


, . . . ( Areal 
Incident wave J 


Area II 


Area III 



(lead 1) 



V = 



v=v„ 



v = o 



FIG. 7: Reflected, refracted and transmitted waves in a beam- 
splitter 

incident wave in lead number one, for example, we can 
calculate the reflected, refracted, and transmitted waves 
in lead 3 and 4, using the continuity of the wavefunction 
and its derivatives at the beamsplitter interface (x=0 and 
x=L). For example, let us start with an incident state 
at energy E with momentum ki in the SO coupling state 
a = +. By the conservation of energy, the reflected wave 
in lead 3 is a superposition of the states \k r (E, +), +) 
and \k r (E, — ), — ) with 



k r (E, +) = h(E, +) = k r (E, -) - Ak. 



(Bl) 



The translational invariance of the beamsplitter along 
the y-axis leads to the conservation of the y-component 
of the momentum 



ki cos 0i — k r (E, +) cos 8 + — k r (E, — ) cos ( 



(B2) 



We deduce that cos6>; = cos6> r+ , but cosf5i cosf? r _. 
There is dispersion due to the SO coupling, leading to 
an angular separation between the + and — states after 
reflection at the beamsplitter (see Fig. ||). This angular 
separation is given by 



r - = arccos[(l 



Afc 



K(E, 



) cos 6i] ^ 9 r 



(B3) 



Starting with an incident state with SO coupling label 
— , the angular separation is 



6 r+ = arccos[(l 



Ak 



k r (E,+) 



) cos 6i] 7^ f5 r _ 



(B4) 



Reflected states 



I ME.-) . - > 




kjCose i+ = k r (E,+) cos G r 
= k (E,-) cos B 



Incident state 



FIG. 8: angular separation after reflection at a beamsplitter 



In analogy with the the total reflection for incident an- 
gles below the critical angle in classical optics, we can 
even have (1 + k ^e+) ) cosOi > 1, leading to a suppres- 
sion of the reflection in the + state for a small enough 
incident angle. We note that starting with a mixture of 
+ and — states in the incident beam of electrons, one 
could suppress the reflection of the + state, thus achiev- 
ing a polarization of the beam. Although interesting, 
we will neglect this effect of angular dispersion by con- 
sidering only non-equilibrium electrons above the Fermi 
energy for which 4p -C 1 (as SO coupling is small com- 
pared to the kinetic energy), and important incident an- 
gles (9i — j). In this case the angular separation is very 
small (Af3 = 6> r _ — 8 r+ < 1). The equations of continu- 
ity of the wavefunction and its derivative are then much 
easier to solve, and one can find that the incident, re- 
fracted, and reflected waves have the same spin at x=0, 
and the refracted and transmitted waves have the same 
spin at x=L. Using Eq.(^) to find the spin overlap be- 
tween different leads we find, starting with a + incident 
state, 



|* r ) = r [cosOi \k,- 



- ism9 t \k, -)] 
(B5) 

where \^t) is the transmitted wave (into lead 4) and 1^,.) 
the reflected wave (into lead 3). As the transmitted wave 
has the same direction of propagation as the incident 
wave, there is no mixing of the SO coupling states, and 
we have only the usual transmission coefficient t. For 
the reflected wave, the direction is changed and we have 
to mix the different SO coupling states to obtain the 
same spin as the incident wave on the interface with the 
beamsplitter. If we start now with a — incident state, we 
find 

|* t ) = t\k,-) |* r ) = r [isinOi \k,+) + cosflj \k, -)]. 

(B6) 

The same analysis can be done for an incident state in 
lead 2 with 9i replaced by — 6^. We then deduce the whole 
scattering matrix: 



S = 



r cos 6 
ir sin t 

t 





ir sin ( 
r cos t 



t 



t 


r cos Qi 
-ir sin f9, 




t 

-irsin6 
r cos 6i 



(B7) 
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